Transition Form Factors and Mixing of Pseudoscalar Mesons from Anomaly Sum Rule 
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I. INTRODUCTION 

The phenomenon of axial anomaly 0, Q plays an 
important role in non-perturbative QCD and hadronic 
physics. The axial anomaly is known to govern the two- 
photon decays of the Tr°,r] and rj' mesons and is usually 
considered for a case of real photons. However, the dis- 
persive form of it Q can be considered for virtual photons 
also leading to a number of interesting applications. 

One of the consequences of the dispersive approach to 
axial anomaly is a so-called anomaly sum rule (ASR) 
It gives, in particular, a complementary way to de- 
scribe the 7T° Q and the rj, vj II, Jl transition form fac- 
tors (later developed also in [10, LUj ) at all Q 2 even be- 
yond the QCD factorization. This is especially impor- 
tant in view of the recent experimental studies of the 
77* — > 7T° (77, 77') transitions [12h14| . In particular, the 
pion transition form factor, measured by the BABAR 
collaboration [l2j], revealed unexpectedly large values in 
the range of Q 2 = 10 — 35 GeV 2 resulting in an ex- 
cess of the pQCD predicted limit 15 1 Q 2 F 7TJ — > \/2f- K , 
JV = 0.1307 GeV. This striking result attracted a lot 
of interest and motivated extensive theoretical investi- 
gations. As a result, the transition form factors were 
(re) investigated using the framework of light cone sum 
rules (l6l - [20) . including the flat- like modifications of the 
distribution amplitude |2lTl23l] , the light-cone holography 
approach [H, [2f| , and various model approaches, like a 
chiral quark model [26} (see also [13413); a vector meson 
dominance model and its modifications [3ll [32j. Some 
other approaches can be found in [33l - l42T |. 

In this paper we extend and develop the anomaly sum 
rule approach 0, H[ to the transition form factors with a 



systematic account of the effects of mixing of rj, rj' mesons 
and quark-hadron duality. Also, we performed a new 
analysis of the anomaly sum rule for 7r° using available 
experimental data, including the most recent ones from 
the BELLE collaboration [lj. 

The paper is organized as follows. In the Section II 
we give an overview of the anomaly sum rule approach 
and apply it to analyze the pion (isovector channel of 
the ASR) and rj, rj' (octet channel of the ASR) transi- 
tion form factors. We paid a special attention to the 
seemingly controversial data from the BABAR [T2J and 
BELLE [l4j collaborations. The analysis of different sets 
of data showed that inclusion of the BABAR data re- 
quires a non- (local) OPE correction to the spectral den- 
sity, while the BELLE data alone neither requires nor 
excludes it. In the Section III we develop and reformu- 
late the description of mixing, which plays a special role 
for the rj, tj 1 mesons, in a way which does not require 
the introduction of intermediate non-physical states. The 
problem of compatibility of different mixing schemes is 
also discussed. In the Section IV we performed a numeri- 
cal analysis of the mixing parameters of the rj — rf system 
basing on the obtained sum rule for the transition form 
factors. The possibility of the non-OPE correction to 
the spectral density in the octet channel is investigated 
as well. The summary is presented in the Section V. 



II. ANOMALY SUM RULE APPROACH 

The triangle graph amplitude, composed of the axial 
current J a § with momentum p = k + q and two vector 
currents with momenta k and q (VVA amplitude) 
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can be decomposed [H} (see also [H, HU) as 
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where the coefficients Fj = Fj(p 2 , k 2 , q 2 ; to 2 ), j = 
1, . . . ,6 are the corresponding Lorentz invariant ampli- 
tudes constrained by current conservation and Bose sym- 
metry. Note, that the latter includes the interchange 
(j, -f-> u, k «-> q in the tensor structures and k 2 f-> q 2 in the 
arguments of the scalar functions Fj. 

In what follows, we consider the case when one of the 
photons is real (k 2 = 0) while the other is real or virtual 
(Q 2 = —q 2 > 0). For the invariant amplitude F 3 — F 6 
the anomaly sum rule (ASR) is given by [j| 
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A ( 3 a) (s,Q 2 - 1 m 2 )ds = ^-N c C^ , (3) 



where A3 — ^Impi{F^ — Fq) 1 N c = 3 is a number of 
colors, to is a quark mass and CW are the charge factors 
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zero in the massless limit [48j|). while the integral 

SZ' 4 0) (s,Q 2 ;m 2 )ds equals exactly ±N C C^ . 

Saturating the l.h.s. of the three-point correlation 
function ([T]) with the resonances and singling out their 
contributions to the ASR © we get the (infinite Q) sum 
of resonances with appropriate quantum numbers 

^fiA = / A a) (s,Q 2 ;m 2 )d S = —N C CV>. 

^ Urn? 27T 

(6) 

Here the projections of the axial currents jj^ onto one- 
meson states M(— 7r°, 77') define the coupling (decay) 
constants 

(0\J^(0)\M(p))=ip a r M , (7) 

while the form factors Fm 7 of the transitions 77* — » M 
are defined by the matrix elements 



of components of the axial currents J^ 5 



d\J kx (M(p)\T{J^x)J v (Q)m = ^v P *k p q°F Ml . 

(8) 

The relation ([6|) expresses the global duality between 
hadrons and quarks. 
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A. Isovector channel (7r°) 

For a case of the isovector channel, the first con- 
tribution is given by 7r°, while the higher contribu- 
tions can be absorbed by the "continuum" contribution 
ffi A 3 3) (s, Q 2 ; m 2 ), so the ASR © takes the form: 
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Our aim is to investigate the transition form factors of 
7T°, 77 and r\' mesons, so hereinafter we restrict ourselves 
to isovector (a = 3) and octet (a = 8) channels of axial 
current. 

For these channels the ASR (O has a remarkable prop- 
erty - both perturbative and nonperturbative corrections 
to the integral are absent [71(. The perturbative cor- 
rections are excluded because of Adler-Bardeen theorem 
[4^ . while the nonperturbative corrections are also ab- 
sent as it is expected from 't Hooft's principle. The 't 
Hooft's principle in its original form [47] implies that the 
anomalies of the fundamental fields are reproduced on 
the hadron level. In the dispersive approach this means 
the absence of the corrections to the dispersive sum rules. 

Let us stress, that the spectral density A^(s, Q 2 ; to 2 ) 
can have both perturbative and nonperturbative cor- 
rections (however, the first-order correction cx a s is 
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*UF^(Q 2 ;m 2 ) = —N C C ( ^ - / Af( S ,Q 2 ;m 2 )ds. 

(9) 

(3) 

The lower limit s of the integral we will refer to as 
a "continuum threshold" , bearing in mind that in a local 
quark-hadron duality hypothesis it has a meaning of the 
interval of duality of pion. Also, it can be determined 
directly from the ASR as it will be demonstrated later. 

The spectral density A3 (s, Q 2 : m 2 ) can be calculated 
from the VVA triangle diagram 
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(10) 

where R{s,m 2 ) = J 1 — ^f-, to is a quark mass. 

From ([9]), (fTU| a straightforward calculation gives an 
expression for the pion transition form factor, 
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where Rq — R(s^\ m 2 ). This expression (to our best 
knowledge, for the first time) takes into account the con- 
tribution of quark mass. However, in what follows, we ne- 
glect the contributions of all light (u, d, s) quark masses 
to the ASR in both isovector and octet channels. The 
phenomenological applications of the quark mass effects 
will be considered in a separate paper. 

In this approximation, the spectral density (|10[) is pro- 
portional to S(s) at Q 2 = 0, so the continuum term in 
the ASR ^ goes to zero. This corresponds to the fact 
that contributions of axial states are zero at Q 2 = 0, 
and contributions of higher pseudoscalar states should 
be suppressed in order for the axial current to conserve 
in the chiral limit. 

Relying on the local quark-hadron duality hypothe- 
sis, the analysis of the two-point correlation function 

(3) 

gives the value for the continuum threshold Sq = 0.75 
GeV 2 [49] . Actually, can be determined directly 
from the high-Q 2 asymptotic of ASR, where the QCD 
factorization predicts the value of transition form factor 
Q 2 F™ = EDI- The high-Q 2 limit of d) (m = 0) 

immediately leads to s 3) = 4tt 2 / 2 = 0.67 GeV 2 . This 
expression, substituted to (fTTj) with m = gives 



2V2tt 2 U 4tt 2 / 2 + Q 2 



(12) 



so it proves the Brodsky-Lepage interpolation formula 
for the pion transition form factor [50| . which was later 
confirmed by Radyushkin 51| in the approach of local 
quark-hadron duality. Let us stress, that in this way we 
found that it is a direct consequence of the anomaly sum 
rule (which is an exact nonperturbativc QCD relation). 

Currently the data from the CELLO [52[, CLEO [Eg, 
BABAR HH and BELLE [3 collaborations cover the 
region of Q 2 = 0.7 - 35 GeV 2 (see Fig. Q). While 
at Q 2 < 10 GeV 2 they are consistent, at larger vir- 
tualities the BABAR and a newly released BELLE 
data are quite different. In this situation we will con- 
sider two sets of data: CELLO+CLEO+BELLE (I) and 
CELLO+CLEO+BABAR (II). 

When compared to the experimental data set (I), Eq. 
(TT21 gives a reasonable, consistent with the data descrip- 
tion (x 2 /d.o.f. — 1.01, d.o.f — 35, see the dashed line in 
Fig. [1]). For the data set (II) the description is worse 
(X 2 /d.o.f. = 2.29, d.o.f. = 37). So, if the data of the 
BABAR collaboration are correct, we come to a violation 
of the ASR-based expression for (|T2"j) . This means 
that the spectral density (|10[) must have a substantial 
correction 5A^, which results in corrections to the con- 



(3) 

J g (3) SA 3 ds and pion 5I„ = J* 4 V, SA 3 ds 
At the same time, as the full integral re- 



tinuum SI cont 
contributions, 
mains constant, the corrections can be related 



It is important, that the main terms of the continuum 
Icont and the pion /„■ contributions have essentially dif- 
ferent Q 2 behaviour 
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so the 1 IQ 2 power correction to the continuum contribu- 
tion is of the order of the main term of the pion contri- 
bution. 

Let us now discuss the sources of possible corrections 
to the spectral densities which in our approach are the 
counterparts of the NLO and higher twist corrections in 
the accurate pQCD fits [l8|. Note that one- loop cor- 
rections to the spectral densities of all structures in the 
VVA correlator in the massless case are zero, which can 
be easily deduced from the results of 1481 . If this nulli- 
fication is due to conformal invariance |54l . |55| , one may 
expect the two- and higher loop corrections to be nonzero 
due to the beta-function effects, which is recently ob- 
served in the soft-photon approximation 56] . Never- 
theless, these higher a s corrections as well as the local 
OPE-induced corrections are small enough to produce 
an enhancement in pion transition form factor shown by 
BABAR: SF^^ ~ log(Q 2 )/Q 2 . Clearly, from dimensional 
arguments, such term cannot appear from the local OPE 
[7j . Thus, to comply with the ASR, the correction should 
be of a non- (local) OPE origin simulating the contribu- 
tion of the operator of dimension 2. Among the possible 
sources of such correction are non-local condensates, in- 
stantons or short strings [571 ]. 

Although the exact form of such correction is not yet 
known, we can construct the simplest form of it relying 
on general requirements. Namely, the correction should 



vanish at s 



(3) 
o 



oo (the continuum contribution van- 
ishes), at s — > (the full integral has no corrections), 
at Q 2 — > oo (the perturbative theory works at large Q 2 ) 
and at Q 2 — > (anomaly perfectly describes pion de- 
cay width). Supposing the correction contains rational 
functions and logarithms of Q 2 , the simplest form of the 
correction satisfying those limits is 



1 a4V n Q 2 

2V2V(4 3) +Q2)2 V s ( 3 ) 
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where A and a are dimensionless parameters. Then the 
pion transition form factor with this correction reads: 



F 7rj (Q 2 ) = ^- (U+6U) 



,(3) 
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For the continuum threshold we will use s 



(3) 



2 f2 



0.67 GeV 2 , which implies that at very high Q 2 the fac- 
torization restores. Nevertheless, if the factorization is 
violated at all Q 2 , one can use a different value for s^. 

(3) 

Moreover, one can even consider the dependence of Sq 
on Q 2 [38(, which can lead to an effective change of a. 

The fit of (JT7D to the data set (II) gives A = 0.14, a = 
-2.36 with x 2 1 d.o.f. = 0.94 (d.o.f. = 35). The plot of 
Q 2 F- K1 for these parameters is depicted in Fig. Q] as a 
solid curve. Note, that the proposed correction changes 
its sign at Q 2 ~ 8 GeV 2 improving the description in 
both regions: at Q 2 < 8 GeV 2 and at Q 2 > 8 GeV 2 . 
Also, F 1X1 (fTT|) with these parameters A, a (obtained from 
the fit to the data (II)) describes well also the data set 
(I) with x 2 /d.o.f. = 0.84 (d.o.f. = 35). 

The summary of the fitting results for F V1 with (Eq. 
(JTrJ) ) and without (Eq. (fT2|)) correction for different sets 
of data are shown in Table I. One can see, that the data 
sets involving the BABAR data require taking into ac- 
count the correction. The data sets which do not involve 
the BABAR data may be described without such correc- 
tion, although the correction may improve description of 
the data. 

Let us emphasize, that the correction (|16p requires a 
log Q 2 term in 5 A3 itself. This form of the correction is 
different from the one proposed in [To| and could match 
it only if the pre-logarithmic factor in (TTB)) did not de- 

(3) 

pend on s . However, such factor would violate the 
mentioned above requirement of nullification of the cor- 
rection (in the limit — > 0). 

Also, there is a clear distinction with the natural emer- 
gence of the log Q 2 term in the triangle amplitude (which 
was used for description of the BABAR data), where tri- 
angle amplitude itself is used as a model for pion transi- 
tion form factor [26j . Such a description is based on the 
consideration of the triangle diagram as a model for the 
transition form factors and applies the PCAC relation for 
the matrix elements involving large virtualities. To our 
opinion this procedure is not justified to the same degree 
of rigor as for the soft processes. 

It is interesting to note, that extending the expression 
for the pion transition form factor (fTTj) into the timclike 

region, one immediately gets a pole at Q 2 = which is 
numerically close to the mass of p meson. This can be a 
kind of interplay between anomaly and vector dominance. 



B. Octet channel (77, r/) 

For a case of the octet channel, where the main contri- 
butions are from rj and 77' mesons, treating the ASR in 
the same way as for the isovector channel leads to 
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1: (color online) . Pion transition form factor (multiplied 
with correction (Eq. (|17p . solid curve) and without 



FIG 

by Q 

correction (Eq. (|12[1 . dashed curve) as a function of Q 2 com- 
pared with experimental data. 



In comparison with the isovector channel, here we have 
some differences. Due to significant mixing in n — n' 
system the 77' meson has a significant contribution to the 
octet channel. Since 77' decays into two real photons, it 
should be taken into account explicitly along with the 77 
meson. 

Let us note, that in the octet channel, the reliable es- 
timation of s from the usual QCD sum rule method 
meets difficulties (see e.g. discussion in Q). Fortunately, 
the ASR approach allows to determine in the octet 
channel from high-Q 2 asymptotic, just the same way as 
in the isovector channel. Generalization of the pion case 
give s the asymptotes for the 77, rj' transition form factors 

MM (M = r,,r]'y. 



Q 2 Ff s 
So, the Q 2 



C^f M 



My =u^~ -j u -JM)- ( 19 ) 

00 limit of the ASR ([T5J gives 
- * (f»,) 2 +2y/2[fif° + ##]). (20) 



In the Q 2 — limit the transition form factors are 
expressed in terms of the two-photon decay widths of 
mesons, so the ASR (TT51) takes the form: 



/«F 7)7 (0) + /«F^ 7 (0) = 



.(8) 



(21) 



where 



4T 



M— >77 



M 



Solving equations (|20|) . ([2T|) with respect to and 
one of the decay constants f^ and substituting them 
into general ASR (TT5)l we can relate the transition form 
factors F r] i 1 (Q 2 ) and the decay constants f%. 



5 





51* = 0: 


x'/d.o.f. 


SI* # : 


x'/d-o.f. A 


a 


CELLO+CLEO+BABAR+BELLE 




1.86 




0.91 0.12 


-2.50 


CELLO+CLEO+BELLE (set I) 




1.01 




0.46 0.07 


-3.03 


CELLO+CLEO+BABAR (set II) 




2.29 




0.94 0.14 


-2.36 


BABAR 




3.61 




0.99 0.20 


-2.39 


BELLE 




0.80 




0.40 0.14 


-2.86 



TABLE I: X 2 /d.o.f. obtained for Eq. JT2) {SI* = 0) and the best fits of Eq. |T7]) {51* / 0) to different data sets. 



Actually, the four decay constants f%j can be related 
basing on a particular mixing scheme. In Q the ASR 
was analyzed for several sets of parameters and mix- 
ing schemes. The more general consideration of mixing 
schemes and extraction of the mixing parameters are per- 
formed in the next sections. 



III. MIXING 

The problem of mixing in the 77 — 77' system is usually 
addressed either in the octet-singlet (SU (3)) or quark- 
flavor mixing scheme (see e.g. [60, [61] and references 
therein). Basically, meson mixing implies that the "non- 
diagonal" decay constants Jm,M — r},rf in Eq. (|Tj) are 
non-zero. The ir° and the isovector current can be decou- 
pled from the 77 — 7/ system and octet and singlet currents 
because of a very small mixing. 

Let us recall the common approach to the mixing, 
when physical states are represented as a linear combi- 
nation of states with definite SU{3)f quantum numbers 
or quark- flavor content. But, since these states do not 
have definite masses, one cannot write the analogue of 
Eq. (O with these states instead of physical states 77, 77'. 
Indeed, if a state has a definite momentum p* 1 it also has 
a definite mass m 2 = p^p^ . 

One can avoid this problem if formulates the mixing 
in terms of the fields 0, related to the physical states \i) 
by the matrix elements (0|(/3j|j) = #y. 

It is well known, that the field of the pion (j>* is defined 
from the divergence of the isovector component of axial 
channel as a PCAC relation, 



(22) 



To consider the 77 and 77' mixing, one can write down 
a straightforward generalization of this relation, 

8^ = FM$, (23) 

where we introduced the matrix notations: 



I/j5 



M 



7 Q . 

J H5 



Ff fa £a fa 
1 •' 1? J if J G 
= I f fP fP 
\Jri J rj' JG 



(<Pv\ 

4>G 



diag{m 2 ,m 2 , ,m%, ■■■). 



(24) 



The vector J p s in the l.h.s. consists of the components 
of the axial current of definite SU{3) symmetry, a so- 
called octet-singlet basis {a — 8, /3 = 0, see Eq. flU), or 
of the components of the axial current with the decoupled 
light and strange quark composition, a so-called quark- 
flavor basis (a = q, /3 = s): 

J% = -^(uinlsu + dj^d), J* s = S7p7 5 s. (25) 

The elements of matrix F are the meson decay constants 
defined in ([7]). The vector <& of physical fields contains 
the fields of 77 and 77' mesons (f> n and (j)^/ and the fields of 
higher mass states, which we denote as <pc. 

It is possible and commonly to get an additional model 
constraint for the matrix F which is fulfilled by apply- 
ing the respective mixing scheme. Let us first introduce 
the new vector of fields 3? relating the SU{3) symme- 
try property or the quark-flavor contents of currents and 
meson fields. The first two components of <& are labelled 
with the same indices a and j3 as the currents and cor- 
respond to the same symmetry {a = 8, /3 = 0) or quark 
content {a = q, /3 = s). The relation between <fr and $ 
is provided by the orthogonal transformation U 



* = u*, # 



0G 



(26) 



V ! / 



In terms of these fields the Eq. (l23l) can be rewritten as 



(27) 



s M j„ 5 = fm$, 



where F = FU, M = U T MU. 

In our notations the octet-singlet (quark-flavor) mixing 
scheme implies that the matrix F has a (rectangular) di- 
agonal form in the respective octet-singlet (quark-flavor) 
basis, 



fa 

fp 



(28) 



This relation can be obtained from the effective la- 
grangian C which contains an interaction term A£.;„ t — 



8 T a 



fa 



sc 



fa^mlkfik, a = a,0. 

k 



(29) 
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Note, that from the requirement that matrix FU has 
a (rectangular) diagonal form (f^5|) immediately follows 
that FF T is a diagonal matrix. So, imposing the mixing 
scheme is equivalent to imposing the constraint for the 
decay constants: 

+ f$f% + ffifa + - = 0. (30) 

Here the sum is over all physical meson states included 
in the vector 

If we restrict ourselves to consideration of the r\ and rf 
mesons only, then the decay constants form a 2x2 matrix 
and in the octet-singlet and quark-flavor bases satisfy the 
respective diagonality constraints, 



cannot be diagonal simultaneously). The bases are com- 
patible only if fs = fo ( f q = f s ), i-e. in case of the exact 
SU(3)f symmetry. 

Although the octet-singlet mixing scheme is more con- 
venient for the anomaly sum rule relations, which are ex- 
act in case of isovector and octet channels, there are argu- 
ments from chiral perturbation theory against it [62l [63[ 
(see also 0,161). 

The mixing scheme-independent extraction of the de- 
cay constants from the experimental data can finally tell 
us which basis is more adequate for describing the mixing 
in the 77 — rf system. This problem will be addressed in 
the next section. 



8 j-0 



J q J 7} * J q' J 'l 
f<? f s _|_ f 1 f s 



0. 

0. 



(31) 
(32) 



For instance, in the case of the octet-singlet mixing 
scheme (I31|) the matrix of decay constants can be ex- 
pressed in terms of one mixing angle 9 and two param- 
eters fs, fo, forming the well-known one-angle mixing 
scheme: 



j8 jr8 

Fso = ( Jb Ih 

J rj J rj' 



fs cos 9 fs sin 9 
-/ sin6> / o cos0 



(33) 



Similarly, if the quark-flavor mixing scheme restriction 
(|32|) is applied, then it is common to express the decay 
constants in terms of parameters cj>, fs, fo, 



( f q f q , 

* qs I « fs 

\ J 7] Jri' 



f q cos <j> f q sin (f) 
-fssmcf) / s cos</> 



(34) 



While either of the mixing schemes (octet-singlet or 
quark-flavor) is self-consistent, they are incompatible [9, 
[SOl ] . Indeed, octet-singlet and quark- flavor bases of axial 
currents are related by means of rotation matrix: 



7 8 
7° 



V(a) 



J H5 



//5 



, V(a) 



/ cos a 
I sin a 



— sma 
cos a 



(35) 

where tana = Then, as follows from ([2^)l . the matri- 
ces of decay constants F Q) g = F (124]) in the octet-singlet 
(f3"3"| and quark-flavor (|3^|) bases are related as 



(36) 



F 80 = V(a)F qs , 



IV. OCTET CHANNEL: NUMERICAL 
ANALYSIS 



In the Section II, as a consequence of the ASR in the 
octet channel, we obtained the relation between transi- 
tion form factors and decay constants of rj and rj' mesons. 
In this section we use this relation to analyze the decay 
constants in different mixing schemes, described in the 
Sec. III. 

First, let us consider the octet-singlet mixing scheme. 
In order to determine the mixing parameters of this 
scheme, we employ the mixing scheme constraint (|31[) 
and the ASR relations (fl8ll. (l20l). (HP). In terms of 



the parameters (|33)) the Eq. (|20[) reads: s, 



47T 2 /,- 



Then, the regions in fs, 9 parameter space which are con- 
strained by the fit of Eq. {T5J {x 2 /d.o.f. < 1) to the 
BABAR data [H| and by the Eq. (J2TJ (the experimental 
errors of r^n-^ are taken into account) are shown in 
Fig. [5J The yellow intersection determines the parame- 
ters, which can be estimated as: fs — (0.88±0.04)/ ff , 9 = 
-(14.2 ±0.7)°. 

In order to determine the constant /q (which does not 
enter the Eqs. (fT8)l . (|20"]) . (f2i~j) in case of the octet-singlet 
mixing scheme) we need an additional constraint. 

As an additional constraint it is convenient to use the 
ratio of radiative decays of J / ^ for which the more cum- 
bersome contribution of the gluonic anomaly is under 
control. Indeed, according to [66j, the radiative decays 
J/ty — » r/(r/')~f are dominated by the non-perturbative 
gluonic matrix elements, and the ratio of the decay rates 
Rj/v = (T(J/*) -> J/*) -> rr/) is given by 



and so 



F 80 Fj = V(a)F qs F q T s V(a) T 



(37) 



We see that in general case the decay constants cannot 
follow the octet-singlet and quark-flavor mixing scheme 
simultaneously (since the matrices FsqFJq an d F qs Fq S 



R 



J/* 



(0 I GG I rf) 



(0\GG\ rf) 



Prf_ 

Pv 



(38) 



where p vW) = M J/9 (1 - m^/Mj^)^. 

The divergencies of the octet and singlet components 
of the axial current in terms of quark fields are 
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FIG. 2: (color online). Octet-singlet mixing scheme parame- 
ters fs,0. Blue region - constraint of Eq. (|18p (y^jd.o.j. < 1), 
red region - constraint of Eq. (|21[> (experimental uncertain- 
ties are taken into account). 




0.6 0.8 1.0 1.2 1.4 



FIG. 3: (color online). Quark-flavour mixing scheme param- 
eters f q , f s - Blue region - constraint of Eq. (1181) (x 2 /d.o.f. < 
1), red region - constraint of Eq. (|21l) (experimental uncer- 
tainties are taken into account). 



<9 M J° 5 = -^(m u n7 5 u + m d d-y 5 d + m s sj 5 s)+ 



1 3a s 
2^ 4tt 



GG. 



(39) 



(40) 



Taking matrix elements of the above expressions between 
vacuum and 77(77') states and neglecting the u and d quark 
masses, the ratio (|38|) can be expressed in terms of 
the decay constants (J7J) as follows: 



The current experimental value of this ratio is Rj/\$ = 
4.67 ±0.15 [___]. 

Employing this ratio for the octet-singlet mixing 
scheme and taking into account Eqs. ([T5]l and ([5T]l one 
can determine the singlet constant: /_ = (0.81 ± 0.07)/,. 
So, the full set of constants of the octet-singlet scheme 
is: 



fa = (0.88 ± 0.04)/,, 
fo = (0.81 ±0.07)/,, 
9 = -(14.2 ±0.7)°. 



(42) 



For the quark-flavour mixing scheme we can perform 
the similar analysis, using the constraint of the scheme 
([32]) . Eqs. ([201), (EH), (gTJ and fitting the ASR (TH) to the 
B AB AR data [13| . The decay constants of quark- flavour 
basis are expressed in terms of those of octet-singlet 



basis / ' , by means of the Eq. (|36|) . In terms of the 
mixing parameters f q ,f s ,<f> (|34| the Eqs. ([20]) . (|4T|) read: 



*< 8 > =(4/3)^(5/* 



R 



(43) 
(44) 



The Eq. (j4"4")l determines the parameter </> = (38.1 ± 
0.5)°. Then the other two parameters f s ,fq can be es- 
timated from the Eqs. jUJ and ([15]). The plot of the 
regions constrained by these equations is shown in Fig. 
131 The red band indicates the constraint of the Eq. (|2"Tj) 
(experimental errors of Rj/<&, ^ ■ n o n ')^"i-y are taken into 
account) and the blue band indicates the fit of the Eq. 
(JTHJ to the BABAR data at ^jd.o.f. < 1 level. One 
can observe two regions where both equations are com- 
patible. We have chosen the physically motivated one, 
where f q ,f s > /,■ So, the yellow intersection in Fig. [__ 
determines the parameters for the quark-flavour mixing 
scheme, which give the following ranges for them: 



/, = (1.20 ±0.15)/,, 
f s = (1.65 ±0.25)/,, 
cf>= (38.1 ±0.5)°. 



(45) 



The obtained mixing parameters (|45[) are i n ag reement 
with those obtained in other approaches [___, [__!, ___} . 

It is interesting also to get the mixing scheme in- 
dependent constraints on the decay constants. Equa- 
tions (f2"Tj) . (|41l) allow to exclude two of the four decay 
constants which enter the ASR (____). The continuum 

In Fig. S] the 
in the space of 



(8) _____ 

threshold s is excluded using (|2"0")l . 
levels of /d.o.f ■ function of Eq. ([1 
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FIG. 4: (color online). Independent (of mixing scheme) esti- 
mation of decay constants. Black thin curve is a \ 2 /d.o.f. = 1 
level, red dot is a minimum of \ 2 of Eq. (|18[) . Green solid 
and orange dashed lines indicate the constraints of the octet- 
singlet (f3lj) and quark-flavor ({32]) mixing schemes respectively. 



constants f%,f2 are shown (BABAR experimental data 
on F lv (~i\(Q ) and mean values of Rjm, r^^/w 2 7 are 
used). One can see, that y^j d.o.f. < 1 requirement 
(black curve) allows a rather wide range of the param- 
eters. However, the minimum of \ 2 is reached at the 
point /8 = 1.11/,,/° = 0.16/, (x 2 /d.o.f. = 0.84, indi- 
cated as a red dot). Then (|2l"j) . (l4"Tj) allow to determine 
the other two constants: = -0.42/,,/°, = 1.04/,. 
Therefore, the full set of decay constants of the mixing 
scheme-independent extraction is: 



Jri J n' 
J rj J rj' 



1.11 -0.42 
0.16 1.04 



(46) 



The constraints of the octet-singlet (green solid curve) 
(J2U and quark-flavor (orange dashed curve) (|3"2")1 are also 
depicted in Fig. 2] We see that both considered mix- 
ing schemes are consistent with the scheme-independent 
analysis based on the ASR at the level of \ 2 1 d.o.f . < 1, 
even if other experimental errors are not taken into ac- 
count. However, the least \ 2 is reached in the region lying 
outside of both mixing scheme curves. Further improve- 
ment of the experimental data can clear up the question 
of validity of different schemes and give more precise val- 
ues of the mixing parameters. 

One can expect that the possible non-OPE correction 
to the spectral density, discussed for the isovector chan- 
nel, should manifest in the same way in the octet channel 
also. Although the BABAR data [13[ do not show such 
a sturdy growth of Q 2 F ril and Q 2 F v >j as they do for 
)J r the octet combination of the transition form fac- 



f®,F v ' 7 ) does reveal a possible growth 



tors Q 2 (f*F v . 
72]. 

Expecting the similarity of the correction to the spec- 
tral density in the isovector and octet channels, we sup- 
pose that the correction in the octet channel has the same 




FIG. 5: (color online). ASR with correction fl|48|). solid line) 
and without correction f (|18[) . dashed line) for the octet-singlet 
mixing scheme parameters (|42[) compared with experimental 
data (blue points - CLEO, red points - BABAR). 



form as (|TB| . 



Sin 



SA^ds 



1 a4 8) q 2 Q 2 

(47) 

This correction results in an additional term SIg, /ir in the 
r.h.s. of ([T8| . 



2^2 4 8 ) + Q2 



+ 



a4V 



2V^VT' 



(4 8) 



+ Q 2 ) 2 



-(In 



Q_ 

J 8 ) 



(48) 



The plot of the Eq. 0SJ (solid line) and Eq. CB 
(dashed line) for different mixing schemes are shown in 
Figs. I5I6I7I The parameters A = 0.14, a = —2.36 are 
taken the same as obtained for the pion case (data set II) , 
the decay constants (central values) for different mixing 
schemes are employed from (|42[) . (|45|) and (|46|) . We see, 
that the current precision of the experimental data on 
the rj and rj' transition form factors could accommodate 
the same correction as in the pion case, however does not 
require it. 

Finally, let us make the following remark. In this 
paper we developed the approach which does not rely 
on introduction of the non-physical states which have 
no definite masses. At the same time, in literature is 
widely discussed a hypothesis (see e.g. [H, [z3|) that the 
transition form factor of the non-physical state \q) = 
-^(\uu) + \dd)) is related to the pion form factor as 

F q ~f(Q 2 ) — (5/3)i 7 ', 7 (Q 2 ) (where the numerical factor 
comes from the quark charges: (e 2 +e 2 )/(e 2 — e 2 ) — 5/3). 
The states \q) and \s) = \ss) are assumed to be expressed 
in terms of the physical states \rj), \r/') via the quark- 
flavor mixing scheme (73|, 



|g) = cos (p\rj) + sin (j)\ri'),\s) = — sin </>| 77)+ sin</>|7/). (49) 
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FIG. 6: (color online). ASR with correction ((|48|). solid line) 
and without correction f (|18|) . dashed line) for the quark- flavor 
mixing scheme parameters (|45|l compared with experimental 
data (blue points - CLEO, red points - BABAR). 



FIG. 8: (color online). Combinations F nl Q 2 (blue solid line) 
and F^i^Q 2 (red dashed line) from Eqs 



tively as functions of Q 



(f5l"j) and (f52j) respec- 
compared with experimental data. 



„ 0.025 
> 



°. 0.020 
"J 0.015 




10 20 30 40 

<£- , GeV 2 



FIG. 7: (color online). ASR with correction ( (J48J) , solid line) 
and without correction ( (|18[) . dashed line) for the mixing 
scheme-independent parameters (146[) compared with exper- 
imental data (blue points - CLEO, red points - BABAR). 



The plot of Eqs. (f5~Tj) . ([52]) with constants from our 
analysis flU f q = 1.20/ w , / s = 1.65/*, <j> = 38.1° in com- 
parison with experimental data is shown in Fig. [8] One 
can observe a reasonably good agreement with the exper- 
imental data. For the decay constants of Ref. [6(| one 
gets a good description also. 

The agreement with the experimental data may indi- 
cate that the effect of strong anomaly for -j^\uu + dd) 
state is small and the strong anomaly predominantly ap- 
pears in the ss channel. This statement can be rigorously 
checked by means of the anomaly sum rule for the singlet 
channel, which we postpone to future work. 



V. SUMMARY 



Then one can relate the form factors: 



—F vl = Frjj cos (j) + Frfy sin c 



(50) 



Let us now try to incorporate this hypothesis into our 
approach. For this purpose, combining (|18l) . ([50]) and 
using (jlip . (|36j) we can express the 77, 77' transition form 
factors in terms of the constants f q , f s , <j>, 



Fr]y (Q ) 



1 s[) 3) (V2/ s cos 4>~ fq si n < 



^ 2 fsh 



4 3) + Q 2 



1 Sq sin ( 



1 (V2fa sin 4> + f q cos < 



(51) 



^ 2 fsU 



4 3) + Q 2 



-1 (8) 

1 Sq COS ( 



^ 2 fs S ^+Q^ 
where S „ 3) = 4tt 2 /. 2 , 4 8) = (4/3)^ 2 (5/, 2 - 2/ s 2 ). 



(52) 



The exact anomaly sum rule allows to derive an ex- 
pression for the pion transition form factor at arbitrary 
Q 2 , giving the proof for the Brodsky-Lepage interpola- 
tion formula. At Q 2 = it is related to the pion de- 
cay width, while at large Q 2 , basing on the factorization 
approach, it allows to define the pion interval of dual- 
ity which numerically appears to be close to the value 
defined from the two-point correlator sum rule analy- 
sis. However, the proposed approach may be applied 
even when the factorization is broken. It was exactly the 
situation supported by BABAR data requiring a small 
nonperturbative correction Q to the continuum spectral 
density. Having dimension 2, it cannot appear in (a lo- 
cal) OPE and should be attributed to, say, instantons or 
short strings. 

In this work we included in our analysis the recent 
BELLE collaboration data. The main conclusions are the 
following. Although the BELLE data themselves may be 
described without the mentioned correction, but they do 
not also exclude its possibility. Unless the BABAR data 
will be disproved, the need for the correction remains. 
This is supported by the Table 1, where the fits for vari- 
ous combinations of the data are shown. 

In the generalization of our approach to the i] and rj 
mesons, the data may be described without such a correc- 
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tion. However, the possibility of the correction, similar 
to the discussed for the pion case is not excluded by the 
current experimental data and supported by the slight 
growth of the octet combination of transition form fac- 
tors. So we can conclude, that the correction to the spec- 
tral density, first introduced in [3] seems to be universal 
for both, isovector and octet channels. This conclusion 
is in an agreement with a recent discussion in fllj ]. 

The special role in the octet channel plays the mix- 
ing. There are two mixing models on the market now: 
the octet-singlet and quark-flavor mixing schemes. We 
reformulate these models without introducing the non- 
physical states with indefinite masses. Each scheme im- 
plies a certain constraint for the meson decay constants 
(|3Tj) . (|32 j) . Both mixing models are compatible only in 
the exact SU(3) / limit. 

Using the data on the transition form factors of the 
r], rf mesons, the ASR allows to extract the set of decay 
constants in the octet-single, quark-flavor schemes as well 
as in the mixing-scheme independent way, if we add an 
additional constraint of ratio of radiation decays of J/ip 
meson (|41l) . It is shown, that the current data precision 



permits both, the octet-singlet and quark-flavor mixing 
schemes. The future improvement of experimental data 
on transition form factors of rj, rj 1 mesons and ratio Rj/^p, 
expected from the BELLE and BES-III collaborations, 
can determine which scheme is more suitable for descrip- 
tion of the mixing in the rj — rj' system. 
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In case of singlet channel (a = 0) this property is violated 
because of gluonic anomaly 

The growth of the octet combination is due to the fact, 
that Q 2 F V1 at Q 2 > 20GeV 2 slightly tends to go up, and 
Q 2 F r]1 tends to go down, while is positive and is 
negative. 

Note, that in our approach the same relation exists only 
for the fields (see Eq. ((26])). 



